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The existing data obtained in deep-inelastic scattering (DIS) experiments in the 
kinematical region xb ~ 10 2 and Q 2 > 5 GeV are examined to gain further insight 
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observed in diffractive scattering and those found in large-transverse-momentum jet 
production, but also the striking features of the structure function F2(xb,Q 2 ) ob- 
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into the dynamics of interacting gluons. It is shown that not only the regularities 



Qh| served in normal DIS-events in this kinematical region can be understood as direct 

i- 

43 



QJ , consequences of self-organized criticality. The results explicitly demonstrate the use- 



fulness of concepts and methods of Complex Sciences in understanding the various 
striking features observed in high-energy collision processes in which "soft" gluons 
play the dominating role. 
Deep-inelastic lepton-nucleon scattering utilizing electron, muon, and neutrino beams 
is an extremely useful tool in studying the structure of matter and its interactions. Since 
the operation of HERA, the first electron-proton collider, the kinematical region of such 
experiments has been significantly extended. In terms of the difference (q) between the 
initial and the final four-momenta (k and k') of the electron, and the four momentum (p) of 
the proton, the kinematical region in which data are now available is: < Q 2 < 9 ■ 10 4 GeV 2 
and 10~ 5 < Xb < 1, where Q 2 = —q 2 , and xb = —q 2 /(2pq)- 

Physics of deep-inelastic lepton-nucleon scattering (DIS) in the "small-x^ region 







(xb ~ 10 -2 , say)" for sufficiently large Q 2 -values plays a very special role. In fact, "Small-xs 
Physics" has become one of the most active fields in Particle Physics since the beginning of 
the 1990's. There are a number of reasons why this field has become so attractive. Some of 
them are listed below: 

(a.l) It has been reported, soon after HERA started to deliver luminosity, by ZEUS- 
CollaborationS and Hl-Collaborationi that, in this kinematical region, the structure function 
F2 rises towards smaller Xb, with the strength of the rise increasing with Q^. Taken together 
with the fact that the valence quarks contribute very little to xb ~ 1CT 2 , the above-mentioned 
experimental resultsEHi imply, not only that there are many gluons in this region, but also 
that the interactions of such gluons play a significant role. Now, since in accordance with 
QED and QCD, only charged particles namely the seaquark(s) and/or antiseaquark(s) can 
be "directly seen" by the virtual photon 7*, the interaction which are responsible for such 
pair-creation and pair-annihilation processes cannot be neglected at the moment when 7* 
"sees something" from the gluon system. This means, in contrast to the kinematical region 
Xb 3> 10~ 2 where preexisting valence quarks dominate, it is not possible to consider as first- 
order-approximation the struck quark or antiquark as a parton (which is by definition free) 
and to improve the approximation by subsequently taking the relevant Feynman diagrams 
into account. This meansi, theoretical assumptions (which specify e.g. which Feynman 
diagrams should be chosen, whether those leading to higher twists should be included, etc.) 
are needed as input — that is before the optimum set of parameters in the gluon PDF's 
(Parton Distribution Function) can be determined by fitting them to the experimental data!. 

(a. 2) It is see ni, that the above-mentioned extension of the parton idea with the help of 
pQCD and empirically determined PDF's indeed work well in reproducing the existing data!. 
Beautiful fits can be found! for all the measurable quantities inside as well as outside the 
small-xs-region (xb ~ 10 -2 ). While such fits must be considered as excellent descriptions of 
the existing experimental data!, questions such as the following can, and should nevertheless 
be raised!: Is this the only way to describe Nature at the basic level of matter where 
experimental information can be obtained by performing high-energy collision processes at 



relatively large momentum transfer (the so-called "hard" or "semi-hard" processes)? With 
such a large number of adjustable parameters in the theory, how much can we learn about 
Nature by comparing data with fitted curves? The reason why such questions deserve to 
be discussed can be readily seen by recalling what has happened in connection with the 
exciting discussions about the possible quark-substructure a few years ago: The excitement 
began when CDF Collaboration observed@ a deviation between their data and a prediction 
based on pQCD and the then best version of PDF's. The excitement ended when a group 
of pQCD- and PDF-experts showed@ that perfect agreement with the data can be achieved 
by readjusting some PDF-parameters concerning the gluons. Having seen this example, we 
can certainly understand why more and more people are now ready to adopt the view-point 
mentioned by Cooper-Sarkar, Devenish and de Roeck in their review article!: "Finally there 
is always the nagging doubt that the freedom to choose a fairly arbitrary function of x for the 
input parton distributions may be hiding the real breakdown of the standard description." It 
is perhaps worthwhile to compare the situation here with a well-known example in a different 
branch of science: While everybody agrees that earthquake is catastrophic, people may 
debate on "What is more useful in understanding Nature? Is it more useful to parameterize 
everthing we know about earthquakes, insert it as input, and make use of the largest and 
fastest computer to "predict" when, where, and how the next earthquake comes? Or, is it 
more useful to discover something like the Gutenberg-Richter law! although it cannot be 
directly used to make "predictions"? 

(/?) Events with large rapidity gaps (LRG) have been observed in this kinematical region!. 
This observation shows that inelastic diffractive scattering can take place not only in hadron- 
hadron collisions, but also in deep-inelastic scattering processes. This observation gives rise 
to a number of questions: Are the mechanisms which lead to such diffractive scattering 
processes in different reactions related to one another? In particular, if it is indeed "the 
exchange of colorless objects" which is responsible for diffractive scattering in hadron-hadron 
collisions, is it the same kind of "colorless objects" which is also responsible for diffractive 
scattering in deep-inelastic lepton-hadron scattering processes!!? What are such "colorless 



objects"? Can the existence of such objects be understood in terms of QCD? What is the 
relationship between such "colorless objects" and the interacting soft gluons mentioned in 



fa)? 



(7) In normal DIS-events, the observed xb- and Q 2 -dependence of F 2 mentioned in (a) 
and the relationship between these two properties can be well-described quantitatively by a 
simple empirical formula proposed by Haidti 

F 2 (x B ,Q 2 ) = a + mlog ^ log — , (1) 
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where the values of the constants are a = 0.074, m = 0.364, xbo = 0.074, Q% = 0.5 GeV 2 . 
As can be readily seen!!, this empirical rule of Haidti indeed gives a very good description 
of the datai in the region xb ~ 10~ 2 , for sufficiently large values of Q 2 (Q 2 > 5 GeV 2 , say)! 
Is this empirical finding merely a parameterization where its simplicity is nothing else but 
a happy coincidence? Should this remarkably simple empirical fact about F 2 (xb,Q 2 ) be 
ignored? Or, is it worthwhile to ask questions such as the following: Can this empirical 
formula be understood in terms of QCD? What is the relationship between this formula and 
the interacting soft gluons which dominate the small-x^ region? What is the relationship 
between this formula and the fact that LRG events exist in the small-xs region? 

The questions mentioned in (a), (/?), and (7) are closely related to one another — so 
close that they should have been discussed at the same time. But, for reasons which will 
be clear later on (see iv, v, and vi below), our research began with the facts and questions 
associated with those discussed in (/3): 

In a recent Letter0 and a subsequent longer paperB, we proposed that the "colorless 
objects" which manifest themselves in LRG events are color-singlet gluonic BTW-avalanches 
due to self-organized criticality (SOC)ffl, and that optical-geometrical concepts and meth- 
ods are useful in examining the space-time properties of such objects. The theoretical 
arguments and experimental facts which support the proposed picture can be summarized 
as follows: 

(i) The characteristic properties of the gluons — especially the local gluon-gluon coupling 



prescribed by the QCD Lagrangian, the confinement, and the non-conservation of gluon 
numbers — strongly suggest that systems of interacting soft gluons are open, dynamical, 
complex systems with many degrees of freedom, and thus in such systems colorless and 
colored gluon-clusters in form of BTW-avalanches (see below) can be formed. 

(ii) It has been observed by Bak, Tang, and Wiesenfeld (BTW)il that a wide class of 
open, dynamical, complex systems, far from equilibrium evolve into self-organized critical 
states, and local perturbations of such critical states may propagate like avalanches caused 
by domino effects over all length scales. Such a long-range correlation effect eventually 
terminates after a time-interval T, having reached a final amount of dissipative energy, and 
having effected a total spatial extension S. The quantity S is called by BTwEl the "size" , 
and the quantity T the "lifetime" of the "avalanche" and/or the "cluster". It is observedEIE 
that there are many such open dynamical complex systems in the macroscopic world, and 
that the distributions D$ of S and the distribution Dt of T of such BTW avalanches/clusters 
obey power laws: Ds(S) oc and Dt(T) oc T~ v , where [i and v are positive real constants. 
Such characteristic behaviors are knowniB as the "fingerprints of SOC". 

(iii) In order to see whether SOC and thus BTW-avalanches can also exist in microscopic 
systems at the level of quarks and gluons, a systematic analysis^! of the data! for diffractive 
deep-inelastic scattering (DIS) has been performed. The results of the analysis can be 
summarized as follows: The "SOC-fingerprints" indeed exist in diffractive DIS in the small- 
xb region (xb ~ 10~ 2 ) where interacting soft gluons play the dominating role. The size and 
the lifetime distributions of the BTW-avalanches observed in such processes indeed show 
power-law behaviors and the exponents are approximately —2. That is 

D S (S) oc S-f, V « 2 , (2) 
D T (T) oc T~ u , v&2. (3) 

which is true in all Lorentz frames. The observed power-law behavior implies in particular 
that colorless gluon clusters are spatiotemporal complexities which have neither a typical 
size, nor a typical lifetime, nor a typical static structure. 



(iv) The usefulness of the proposed SOC picture has been demonstrated in |TT] and in 
||12|| . It is shown in particular that simple analytical formulae for the differential cross- 
sections da/dt and d 2 a/dtd(M%/s) can be derived for inelastic diffractive scattering, not 
only for small-x^ DIS and for photo-production but also for proton-proton and antiproton- 
proton collisions. It has been pointed out that color- singlet gluon clusters (cq) can be 
readily examined00 experimentally in inelastic diffractive scattering processes, because the 
interactions between the struck Cq and any other color- singlets are of Van der Waal's type 
which are much weaker than color forces at distances of hadron-radius. Thus not very much 
momentum need to be transfered to a Cq by the projectile (which can be a 7*, a 7, a p 
or a p) in order to "knock it out of the mother proton". For this reason, by considering 
inelastic diffractive scattering0'@, we were able to check the existence and the properties of 
the color-singlet gluon-clusters. 

(v) After having seen0ll! the existence of SOC-fingerprints, and having demonstrated0ll! 
the usefulness of such concepts and methods in Particle Physics by confronting them with 
the large body of experimental data@ on inelastic diffractive scattering where the color- 
singlet and only the color-singlet and only the color-singlet gluon-clusters can be calculated, 
we discussed as the next step, the question: "What can we say about those gluon-clusters 
which are not color-singlets?" Here, we have to recall that, due to the (experimentally 
observed) SU(3) color-symmetry, most of such gluon-clusters are expected to carry color 
quantum numbers; that is, they are color-multiplets (instead of singlets) which will hereafter 
be denoted by c*'s. On the other hand, in accordance with the (experimentally confirmed) 
characteristic features of the BTW theoryEl'0, the existence of SOC fingerprints in gluon 
systems cannot, and should not, depend on the dynamical details of their interactions - 
in particular not on the intrinsic quantum numbers they exchange during the formation 
process. This means in particular, that the size (S) distribution Ds{S) and the lifetime (T) 
distribution D T (T) of the colored gluon clusters (the c*'s) should not only exhibit power- 
law behavior, but also have the same power as that found for color-singlet counterparts 
in diffractive scattering processes0Jii! Is the existence of SOC in open dynamical complex 



systems indeed so general, and its characteristic features indeed so universaP. Can such 
expectations be checked experimentally! 

(vi) The questions raised in item (v) has been discussed in a more recent LetterS where in 
particular the following has been pointed out: Similar to individual quarks (g's) or antiquarks 
(g's), colored gluon-clusters (c*'s) can also be "knocked out" of the mother proton p by a 
projectile, provided that the corresponding transfer of momenta is large enough. However, 
in contrast to the knocked-out g's or g's (in usual DIS-events), the knocked-out c*'s may 
or may not have "color lines" connected to the remnant of the proton — depending on 
the color-quantum number carried by the final state of the struck c*, after it absorbs the 
projectile which in proton-proton (pp) or proton-antiproton (pp) collisions can be a quark 
or an antiquark or a gluon. Such knocked-out c*'s manifest themselves in form of hadronic 
jets. In fact, it has been explicitly shown in@ that the observed high-i^y-jets in pp collisions 
can be described by the proposed SOC-picture. The same idea and the same method can 
be applied to lepton-induced DIS in which the large momentum transfer is delivered by the 
virtual photon 7*. The obtained results will be shown elsewhere. 

We note in deriving the results presented mS the following properties of the colored 
gluon clusters in form of BTW- avalanches have been explicitly taken into account: 

First, because of the universality and the robustness of SOC, the formation processes and 
the properties of the BTW-avalanches, in particular the SOC-fingerprints are expected to 
be independent of the intrinsic quantum number they carry. Having these and the following 
arguments in mind, we conclude that the size- and the lifetime-distributions of the colored 
gluonic BTW-avalanches are expected to be the same as those for colorless ones which have 
already been experimentally examined in diffractive scattering. It is useful to note that, 
inside a BTW-avalanche every constituent is in general interacting with more than one of 
its neighboring constituents, in fact, everyone of them is interacting directly or indirectly 
with everyone else through color forces. Hence, interaction between a BTW-avalanche and 
an incident space-like photon 7* (or a hadron h) is weaker than the (average) interaction 
between the constituents of the BTW-avalanche. This means, in a collision process between 
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7* (or h) and a BTW-avalanche the latter acts as entire object — independent of the fact 
whether the interaction between the struck BTW-avalanche c* and its neighbors is of Van 
der Waal's type. In other words, the question whether a BTW-avalanche (either colorless 
or colored) can be knocked-out depends solely on the momentum-transfer it obtains in the 
collision process. 

Second, since gluonic BTW- avalanches are spatiotemporal complexities which have nei- 
ther a typical size, nor a typical lifetime, nor a static structure, the scattering process 
between the virtual photon 7* and the struck avalanche c* is not a process in which "the 
electromagnetic structure of that avalanche" is being probed! The 7* interacts within the 
allowed time interval with more than one of the charged constituents of c*, where the dynam- 
ical details about the individual subprocesses are rather unimportant. Roughly speaking, 
the role played by 7* is simply to deliver a sufficient amount of momentum-transfer to the 
gluonic avalanche c* such that the entire c* can be knocked-out of the mother proton. 

Let us now turn our attention to the majority of normal DIS-events, in which the virtual 
photon 7* encounters a gluonic avalanche c*. We note, based on the facts mentioned above, 
it is important to know the probability for a 7* to meet such a c* of size S and lifetime T; and 
this can be immediately written down as the product of Ds{S)Dt{T) and ST which is the 
space-time volume of the spatiotemporal complexity c*. Since F 2 {xb ~ 10~ 2 , Q 2 > 5 GeV 2 ) 
is in fact nothing else but the total probability for the above-mentioned interaction in the 
given kinematical range to take place, we need to collect all those terms which contribute 
to the total probability. In doing so, we are led to the conclusion: 



where D${S) and Dt{T) are given by Eqs.(^) and (|3|) respectively, N is a normalization 



F 2 (x B ,Q 2 ) 




(4) 



constant, and the integration limits are functions of xb, Q 2 , and P =\P 



max 



xboP, 



'mm 



X B P, 



(5) 
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4P x B rp _ 4P x B 
Ql 1 - •'•/-' """" Q 2 1 ./•/>• 



2 i ' -^min ^15, 1 • (6) 



The yet undetermined constants xbo and Qo can be estimated theoretically (see below). 

The facts and arguments which lead us to Eqs.@ and (|j) are the following: Since 
whatever the charged objects hit by 7*(£_b, Q 2 ) may be, they are part of one of the BTW- 
avalanches which dominate the small-x^ region. S max and S m i n is proportional to the max- 
imum and the minimum amount of dissipative energy this particular BTW-avalanche can 
carry. That is, they are given by Eq.(||), together with xbo ~ 10~ 2 Jl7|]. From Eqs.([4[) 
and (|^) we see, as expected, that the largest contribution for F 2 (xb,Q 2 ) comes from the 
avalanches with the smallest size. Next, we consider the interaction-time r int of such a 
collision event in a proper (e.g. cm.) Lorentz-frame: This can be estimated with the 
help of the Uncertainty Principle by calculating q° (of q = k — k'), where the result idHlll! 
Ti n t = = 4|P|Q _2 xb(1 — xb)~ x ■ The lower limit of T is determined by the requirement 
that the encountered BTW-avalanche has to live long enough to been "seen" by the virtual 
photon ^(xb, Q 2 )- The upper limit T max is determined by the requirement that the charged 
object (s) which carries (carry) xb has to be recognizable by 7* (in the sense that 7* should 
be able to find out whether they are part of the gluon-cluster of size S and lifetime T). 
Hence the resolution power (in the transverse directions, 1/Q 2 ) of 7* has to be sufficiently 



large (that is Q 2 > Q 



By inserting Eqs.(Q), (|3|), (|5|), and (g) into Eq.(|J) we obtain Eq.(|I]) where m is related 
to the normalization constant iV in Eq.([|), and a can be interpreted as the averaged value 
of the contributions from the valence quarks in the small-Xs regionlil. In other words, 
the empirical formula proposed by Haidti, as shown in Eq. ([!]), turns out to be a natural 
consequence of the proposed0S0 SOC-picture for interacting soft gluons which dominate 
small-x^ deep- inelastic scattering processes. However, while Haidt, from an experimental 
point of view, tried to extend the formula ([!]) to (if possible) all values of xb, we see here, 
that, from a theoretical standpoint, it is expected to be valid only in the small-x^ region. 
The comparison with the experimental data is shown in Fig. |1|. 
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Fig. 1. The data for the proton structure function Fz{xb, Q 2 ) are taken from [10|. The solid lines 
are the calculated results by using the Haidt-formula H shown in Eq.(|]). According to the theory 
presented in this paper, the formula Eq.(ffl) should only be valid for xb~ 10~ 2 . This is indicated 
by vertical dashed lines. The Q 2 values corresponding to each box are given in brackets. 
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